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ABSTRACT 

GRO J2058+42 is a 195 s transient X-ray pulsar discovered in 1995 with 
BATSE. In 1996, RXTE located GRO J2058+42 to a 90% confidence error circle 
with a 4' radius. On 20 February 2004, the region including the error circle was 
observed with Chandra ACIS-I. No X-ray sources were detected within the error 
circle, however, 2 faint sources were detected in the ACIS-I field-of-view. We 
obtained optical observations of the brightest object, CXOU J205847. 5+414637, 
that had about 64 X-ray counts and was just 0.3' outside the error circle. The op- 
tical spectrum contained a strong Ha line and corresponds to an infrared object 
in the 2MASS catalog, indicating a Be/X-ray binary system. Pulsations were 
not detected in the Chandra observations, but similar flux variations and dis- 
tance estimates suggest that CXOU J205847.5+414637 and GRO J2058+42 are 
the same object. We present results from the Chandra observation, optical ob- 
servations, new and previously unreported RXTE observations, and a reanalysis 
of a ROSAT observation. 

1 

Subject headings: accretion — stars:pulsars:individual:(GRO J2058+42) — X-rays: 
binaries 



1. Introduction 

1.1. Be/X-ray binaries 

The majority of the known accretion- 
powered pulsars are transients in binary 
systems with Be (or Oe) stars (Appa- 
rao 1994; Coe 2000). Be stars are main- 
sequence B stars showing emission in the 
Balmer lines (Porter & Rivinius 2003; Slet- 
tebak 1988). This line emission, and a 
strong infrared excess in comparison to 
normal stars of the same stellar type, is as- 
sociated with circumstellar material which 
is being shed by the star into its equatorial 
plane. The exact nature of the mass loss 
process is unknown, but it is thought to 
be related to the rapid rotation, which is 
typically near or above 70% of the critical 
break-up velocity (Porter 1996; Townsend, 
Owocki, & Howarth 2004). Near the Be 
star the equatorial material probably forms 
a quasi-Keplerian disk (e.g., Quirrenbach 
et al. 1997; Hanuschik 1996), believed to 
fuel the X-ray outbursts. 

In the "normal" outbursting behav- 
ior (type I) of a Be/X-ray binary, series 
of moderate luminosity (10 35-37 erg s _1 ) 
outbursts occur, with each outburst near 
the periastron passage phase of the sys- 
tem's wide eccentric orbit (Stella et al. 
1986; Bildsten et al. 1997). Some systems 
also show infrequent "giant" (type II) out- 
bursts, with luminosities of 10 38 erg s _1 or 
greater. In most systems there is no clear 
correlation between X-ray outbursts and 
optical activity within single outbursts. 
The X-ray activity however, follows the 
long-term optical activity cycle of the Be 
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star, in the sense that no outbursts occur 
in periods where optical indicators of the 
Be star disk, such as Ha emission, have 
disappeared. Periods of X-ray quiescence 
when the Be disk is present are also ob- 
served; these may be due to the truncation 
of the Be star disk well within the neutron 
star orbit (Negueruela & Okazaki 2001; 
Negueruela et al. 2001). 

1.2. GRO J2058+42 

GRO J2058+42 was discovered as a 
198 s pulsar with the Burst and Tran- 
sient Source Experiment (BATSE) on 
the Compton Gamma Ray Observatory in 
1995, during a giant outburst (Wilson et 
al. 1995, 1998). In the giant outburst, the 
pulsar spun-up from a period of 198 s to 
a period of about 196 s. The giant out- 
burst was followed by a series of normal 
outbursts, spaced at 55-day intervals. For 
the first 11 normal outbursts, there was 
an alternating pattern in 20-50 keV pulsed 
flux (i.e. odd outbursts, counting from the 
giant outburst, were brighter than even 
outbursts). However, the All-Sky Moni- 
tor (ASM) on the Rossi X-ray Timing Ex- 
plorer (RXTE), which observed the phase 
averaged flux in the 2-10 keV band for 10 
of these outbursts, did not see such a pat- 
tern, suggesting that the pattern may have 
been caused by changes in the spectrum or 
pulsed fraction. Later, the outbursts con- 
tinued, but the odd-even pattern stopped 
in the 20-50 keV pulsed flux. All of the out- 
bursts showed frequency increases which, 
if attributed to mainly orbital effects, re- 
quired an orbital period of 55 days (Wilson 
et al. 1998; Wilson, Finger, & Scott 2000). 
The outburst behavior, specifically the gi- 
ant outburst followed by a series of normal 
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outbursts, suggested a Be star companion. 

Initially, GRO J2058+42 was approxi- 
mately located with BATSE (Wilson et al. 
1995) and with the Oriented Scintillation 
Spectroscopy Experiment (Grove 1995). 
In 1996 November, scanning observations 
with the RXTE Proportional Counter Ar- 
ray (PCA) further reduced the error region 
to a 90% confidence 4' radius error circle 
centered on a = 314.75°, 5 = 41.72° (Wil- 
son, Strohmayer, & Chakrabarty 1996). 
Castro-Tirado & Birkle (1996) reported 
two optical objects, with magnitudes 18.5 
and 19.5, respectively, in the RXTE error 
circle at a = 314.7417°, 5 = +41.7183° 
and a = 314.7821°, 5 = +41.7461° (J2000, 
uncertainty ±1") on October 15, 1996. 
However, we found no X-ray object at 
either position with ROSAT or Chan- 
dra . The ROSAT HRI catalog, avail- 
able through HEASARC (ROSATHRITO- 
TAL) 1 and through MPE (1RXH cata- 
log) 2 , lists 7 sources in the ROSAT HRI 
field of view from a 1997 June 23 observa- 
tion centered on the RXTE position, but 
only 3 of these sources were unique and 
none were within the RXTE error circle. 
Recently, Reig, Kougentakis, & Papamas- 
torakis (2004) performed an optical pho- 
tometric and spectroscopic analysis of the 
field around GRO J2058+42 from the Ski- 
nakas observatory and suggested the star 
approximately located at a — 314.697°, 
5 = +41.777° (J2000) as its optical coun- 
terpart. 

In this paper we report on Chandra ob- 
servations, archival and new RXTE obser- 



1 http://heasarc. gsfc.nasa.gov/W3Browse/rosat/ 
roshritotal.html 

2 http: / /www. xray.mpe.mpg.de /cgi-bin/rosat / src- 
browser 



vations, archival ROSAT observations, and 
give the details of the optical observations 
and analysis. The combination of archival 
and new X-ray and optical data allowed 
us to pin down the X-ray source and iden- 
tify its likely optical counterpart with a 
moderately reddened Be star. We discuss 
the possible association between a Chandra 
source and GRO J2058+42. 

2. GRO J2058+42 X-ray Observa- 
tions with RXTE 
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Fig. 1.— RXTE ASM data for GRO 
J2058+42 divided in 4 equal ~ 2.1 year in- 
tervals (1996 February-1998 March, 1998 
March-2000 May, 2000 May-2002 June, 
and 2002 June-2004 July) and folded at the 
55.03 day period reported in Wilson, Fin- 
ger, & Scott (2000). Filled circles denote 
> 3(T detections. 

Figure 1 shows all-sky monitor (ASM) 
data from the entire RXTE mission through 
2004 July divided into 4 equal ~ 2.1 year 
intervals and epoch-folded using the em- 
pheris T = MJD 50411.3 + 55.03iV (Wil- 
son, Finger, & Scott 2000). From 1996- 
2002, GRO J2058+42 was detectable with 
the ASM from about 3-5 days before until 
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about 6-8 days after the predicted outburst 
peak. 

In early 1998, two outbursts of GRO 
J2058+42 were observed with the RXTE 
PCA and High Energy Timing Experiment 
(HEXTE). Fourteen observations were ob- 
tained from 1998 January 23 - February 
24 and 13 observations from 1998 March 
18-30. Details including start and stop 
times for each observation can be found 
through the HEASARC 3 . These observa- 
tions spanned from about 16 days until 2 
days before the outburst peak. Pulsations 
were detected from about 6.5 days before 
the peak until the observations ceased with 
r.m.s. pulsed fluxes of 0.2-1.4 mcrab over 
the full 2-60 keV PCA energy range. 

We analyzed energy spectra from archival 
RXTE data for these two outbursts to look 
for spectral variations with intensity. Pre- 
vious efforts (Wilson-Hodge 1999; Wilson, 
Finger, & Scott 2000) reported only on 
spectral analysis of the brightest of these 
observations. We analyzed archival PCA 
Standard2 data and HEXTE science event 
mode E_8us_256_DXlF data for all 27 of 
the observations, using FTOOLs 4 version 
5.3 (Blackburn 1995) to generate spec- 
tra and response files. Background spec- 
tra for the PCA observations were gener- 
ated using the "faint" source models. In 
XSPEC 5 version 11.3.1 (Arnaud 1996), the 
fits included PCA data from 2.7-25 keV 
and HEXTE data from 11-50 keV. All of 
the observations were well fitted with an 
absorbed thermal Bremsstrahlung model 

3 http: / /heasarc. gsfc.nasa.gov 
4 http: / /heasarc. gsfc.nasa.gov/ftools 
5 http://heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
index.html 



(phabs 6 *bremss) with reduced x 2 values 
in the range 0.8-1.3 for 126 degrees of 
freedom. More sophisticated models, such 
as a power-law with a high-energy cutoff, 
were not used because the Bremsstrahlung 
model adequately fit the data and pro- 
vided a simple measure of spectral hard- 
ness across all source intensities. The data 
were better fitted with the Bremsstrahlung 
model than with a power-law at higher 
intensities. Estimated 2-50 keV fluxes 
ranged from 6.8 x 10~ 13 erg cm -2 s" 1 to 
2.6 x 10~ 10 erg cm -2 s _1 . From about 8-9 
days before the predicted peak, the flux be- 
gan to increase, rising above a few xl0~ 12 
erg cur 2 s" 1 . The spectra were well de- 
termined for fluxes above about 2 x 10~ n 
erg cm -2 s -1 , from about 6.5 days before 
the predicted peak, when pulsations were 
also detected. The absorption, Ah, ap- 
peared to be constant with best fit values 
in the range (4.6 - 5.4) x 10 22 cm" 2 . The 
temperature, kT, increased as the inten- 
sity increased, with best fit values from 
10.3±0.5 keV at 2.9 x 10~ n erg cm" 2 s" 1 to 
22.2 ±0.4 keV at 2.6 x 10" 10 erg cm" 2 s" 1 . 

Because of the 1° FWHM field-of-view 
of the PCA, these spectra included emis- 
sion from both GRO J2058+42 and the 
Galactic ridge. The diffuse emission from 
the Galactic ridge is not well described 
at high Galactic longitudes (for GRO 
J2058+42 £ = 83.6°). In fact, the models 
of Valinia & Marshall (1998) and Revnivt- 
sev (2003) predict a Galactic ridge flux of 

6 Abundances are given by Anders & Grevesse 
(1989) and cross-sections by Balucinska-Church & 
McCammon (1992) with a new He cross-section 
based on Yan et al. (1998) in the phabs model, 
used in fits to RXTE and Chandra data. Thomp- 
son scattering is not included in the phabs model. 
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~ 2 x 10~ n erg cm" 2 s" 1 (2-10 keV), tak- 
ing into account GRO J2058+42's galactic 
latitude (b = —2.6°), an order of magni- 
tude higher than our faintest flux mea- 
surements for the GRO J2058+42 region. 
The column density measured with RXTE 
was nearly an order of magnitude larger 
than the Galactic value of 7 x 10 21 cm -2 
calculated with Colden, the Galactic Neu- 
tral Hydrogen Density Calculator 7 , based 
on Dickey & Lockman (1990), suggesting 
there was a significant column of material 
intrinsic to the GRO J2058+42 system. 

In preparation for our Chandra observa- 
tions, RXTE observed the GRO J2058+42 
region in a series of 27 observations from 
2003 December 18 - 2004 January 1, span- 
ning from 6.5 days before until 7.5 days 
after the predicted peak and in a second 
series of five very short (< 1 ks) observa- 
tions from 2004 February 15-19, from 2.5 
days before until 1.5 days after the pre- 
dicted peak. Pulsations were not detected 
in any of the 2003-2004 observations at or 
above the previous minimum detection of 
an r.m.s. pulsed flux of 0.2 mcrab in the 
full 2-60 keV band. Most of the observa- 
tions in 2003 December were short, < 3 ks. 
For the eight longer observations from 2003 
December, with exposure times of 6-14 ks, 
we fitted 2.7-15 keV PCA data in XSPEC 
with an absorbed thermal Bremsstrahlung 
model, with N H fixed at 2.1 x 10 22 cm" 2 , 
the best-fit value measured with Chandra 
(See Section 3.1.2). Best fit temperatures 
ranged from 3-5 keV and 2-10 keV fluxes 
from (2 — 3) x 10~ 12 erg cm~ 2 s -1 , with 
no evidence of a correlation. These fluxes 
should be treated as upper limits for the 
GRO J2058+42 flux, since no pulsations 

7 http: / / asc.harvard.edu/toolkit /colden.jsp 



were detected with RXTE and the spectra 
also include Galactic ridge emission. 

3. X-ray Observations of 

CXO J205847.51+4146373 

3.1. Chandra 

On 20 February, 2004, Chandra ob- 
served the field containing GRO J2058+42 
for 9.9 ks with the Advanced CCD Imaging 
Spectrometer (ACIS) array and the High 
Energy Transmission Grating (HETG) in 
the faint, timed-exposure mode with a 
frame time of 3.241 seconds. The Imaging 
Array (ACIS-I) field-of-view was centered 
on a = 314.75°, 5 = +41.72° (J2000), the 
center of the RXTE error circle for GRO 
J2058+42. The ACIS-I chips were on (10, 
II, 12, 13) and chips S2 and S3 were also 
active. Standard processing v. 7. 7.1 ap- 
plied aspect corrections and compensated 
for spacecraft dither. Level 2 events were 
used in our analysis. Data in the energy 
range 0.5-8 keV were used for all analyses 
to reduce background. We simultaneously 
processed these data using Chandra Inter- 
active Analysis of Observations (CIAO) 8 
version 3.1 and CALDB 9 version 2.27 and 
with the techniques described in Swartz et 
al. (2003). 

3.1.1. Source Detection 

We used the same source finding tech- 
niques as described in Swartz et al. (2003) 
with the circular- Gaussian approximation 
to the point spread function, and a mini- 
mum signal-to-noise (S/N) ratio of 2.6 re- 

8 http : / /cxc . harvard . edu/ciao 

9 See http://asc.harvard.edu/CIAO/ for more infor- 
mation 
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Table 1 

Sources Detected with Chandra 



NAME 


R.A. 


Dec 


7 1 


N b 


S/N c 


4 


USNO e 


2MASS 6 


ROSAT e 




J2000 


J2000 


(") 




(") 


(") 


(") 


(") 


101 


314.67230° 


41.58780° 


10.6 


10 


2.9 


4.2 








CXOU J205847.5+414637 


314.69794° 


41.77704° 


3.7 


64 


7.3 


1.4 


0.37 


0.35 


0.48 



a Extraction radius, defined as 2.5<r, where a is defined in Section 3.1.1. 
b Approximatc number of source counts. 
c Detection signal-to-noise ratio. 

d X-ray position uncertainty (99% confidence radius). 

c Radial separation between X-ray position and cataloged position of counterpart. 



suiting in much less than 1 accidental de- 
tection in the field. The corresponding 
background-subtracted point source detec- 
tion limit is ~10 counts which corresponds 
to a flux (0.5-8.0 keV) of about 7 x 10~ 15 
erg cm -2 s _1 for an unabsorbed power-law 
of spectral index —1.5. Two sources were 
detected. Table 1 gives the source posi- 
tions, the extraction radius, the net counts, 
the signal-to-noise ratio, and the associ- 
ated uncertainty in the X-ray position. Ta- 
ble 1 also indicates whether or not the 
source is identified with a counterpart in ei- 
ther the United States Naval Observatory 
Catalog (USNO-B1.0, Monet et al. 2003), 
the Two Micron All Sky Survey (2MASS) 
and the various ROSAT catalogs and lists 
the angular separation between the X-ray 
position and the optical/infrared/iZOSVlT 
positions. 

The 99% confidence positional uncer- 
tainty listed in column 7 of Table 1 is given 
by r = 3.03(cr 2 /AT + a 2 Q + al) 1 / 2 where 
a is the standard deviation of the circu- 
lar Gaussian that approximately matches 
the PSF at the source location, N is 
the aperture-corrected number of source 



counts for a 2.5 a extraction radius, a 
represents the uncertainty in the image re- 
construction produced by irremovable con- 
tributions to the relative aspect solution, 
and 0\ accounts for the uncertainties in 
the absolute aspect solution. These latter 
two uncertainties are discussed in Weis- 
skopf et al. (2003) 10 . We have used 0.2"as 
a conservative estimate for a Q , and 0.38" 
for <j\. 

The standard CIAO tool wavdetect was 
also used to search for sources using an 
input event file binned by a factor of 4. 
The 2 sources found with wavdetect were 
consistent with the sources listed in Ta- 
ble 1. Both sources are outside the 90% 
confidence 4' RXTE error circle for GRO 
J2058+42. CXOU J205847.5+414637 is 
just 18" outside the error circle, while 101 
is 4.5' from the error circle. 



See also http://cxc.harvard.edu/cal/ASPECT/celmon/ 
and http://cxc.harvard.edu/cal/ASPECT/img_recon/ 
report.html 
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3.1.2. X-ray Spectral Analysis 

We show the X-ray "color-color" dia- 
gram in Figure 2. Since both sources 
were detected with only a small number of 
counts, the uncertainties are large and it is 
difficult to draw firm conclusions, although 
it is clear that CXOU J205847.5+414637 
is very hard and most of the detected flux 
appears above the iridium m-edge (circa 2 
keV) in the telescope response. 

To extract a spectrum and associated 
response files for CXOU J205847.5+414637, 
we used a source region consisting of a 
5.5" radius circle centered on the best- 
fit CXOU J205847.5+414637 position. A 
5.5" circle corresponds to the radius that 
encloses 99.9% of the circular Gaussian 
that approximately matches the PSF at 
the source location. (See Section 3.1.1.) 
We extracted our spectrum from a fits file 
already filtered to contain only events with 
energies between 0.5-8 keV, with the CIAO 
tool psextract. The small number of pho- 
tons detected did not allow us to use the 
X 2 statistic, so we used the Cash statistic 
(Cash 1979) in XSPEC 11.3.1. Since Cash 
statistics cannot be used on background 
subtracted spectra, we did not subtract 
a background; however, the errors intro- 
duced from this are likely to be small be- 
cause we expect about 0.7 background pho- 
tons in our source extraction region. We 
fitted an absorbed power-law spectrum to 
the 64 Chandra counts. The resulting fit 
parameters were: Nr = (2.1 ± 1.0) x 10 22 
cm" 2 and photon index = 1.8 ±0.6. Using 
the error option on the flux command in 
XSPEC, we estimated that the 90% confi- 
dence range for the 2-10 keV unabsorbed 
flux was (3 — 9) x 10~ 13 erg cm" 2 s _1 . The 
quality of the fit was investigated by gen- 




-0.5 0.5 1 

(H-M)/(H+M+S) 



Fig. 2. — X-ray "color-color" diagram for 
the two sources detected with the ACIS- 
I array. The energy bands S, M, and H 
are as follows: S = (0.5 - 1.0) keV, M = 
(1.0 - 2.0) keV, and H = (2.0 - 8.0) keV. 
The solid lines are contours for power-law 
spectra of constant photon number index 
ranging from -1 (innermost) to -4 (outer- 
most) where N H is varying. The dashed 
lines are contours of constant N H for a 
power-law spectrum of varying spectral in- 
dex. Nfj is 0.1, 1, 2, and 5 xlO 21 cm" 2 
from the innermost to the outermost con- 
tour. Thus a source with spectral index -1 
and N H of 10 20 cm -2 would be placed on 
the plot at the intersection of the dashed 
and solid lines at approximately (0.0, 0.1). 
The very hard source to the right is CXOU 
J205847.5+414637 and the source farther 
left is IOl. 

erating 10,000 Monte Carlo simulations of 
the best-fit spectrum. If the model fit the 
data, the fit quality was ~ 50%, mean- 
ing that approximately 50% of the simula- 
tions had Cash statistics lower than that 
of our data. For the power-law model, the 
fit quality was 60%. We also fitted the 
data with a thermal models, a blackbody 
(phabs*bbodyrad) and a Bremsstrahlung 



7 



(phabs*bremss), but in both cases the fit 
quality was poorer, about 70% and 80%, 
respectively. For the blackbody, the emit- 
ting region was very small, with a radius 
of ~ 0.1 km. 

3.1.3. Timing 

We searched for, but did not detect, 
pulses from GRO J 1958+42 using events 
extracted as discussed in Section 3.1.2. 
The search was conducted in the frequency 
range of 5.045-5.155 mHz using the Z\ test, 
with the number of Fourier terms, n, rang- 
ing from 1 to 6. The maximum values of 
found were all consistent with no pulses be- 
ing present. Monte-Carlo simulations were 
used to determine 95% confidence upper- 
limits to the pulse fraction of 37% rms for a 
simple sinusoidal pulse, and 43% rms for a 
pulse profile with Fourier power limited to 
the pulse frequency and its first harmonic. 
For comparison, during the outbursts de- 
tected with the RXTE PCA, the 2-10 keV 
pulsed fractions were 15-24% rms. 

3.2. Archival ROSAT Observations 

Using a PSF-fitting algorithm, we ana- 
lyzed archival data from a 2.4 ks ROSAT 
High Resolution Imager (HRI) observa- 
tion, centered on the RXTE position for 
GRO J2058+42, performed on June 23, 
1997. Based on the outburst ephemeris 
of Wilson, Finger, & Scott (2000), this 
observation took place about 8.5 days be- 
fore the predicted GRO J2058+42 out- 
burst peak. No bright transient object 
was detected. A total of 5 X-ray sources 
were detected in the full HRI-FOV and are 
listed in Table 2. None of these was in- 
side the 4' error circle. The object near- 
est to the error circle was object # 2, and 



was 0.4' outside the error circle. This, as 
well as the other 4 sources, had count rates 
which were below the sensitivity thresh- 
old of the ROSAT all-sky survey (RASS) 
done between June 1990 to January 1991 
(when no source was detected). Objects 
#1,2, and 5 were also found with the stan- 
dard ROSAT analysis and are listed in the 
1RXH catalog. However, the two sources 
consistent with object #1 were flagged as 
close to the detector structure and possi- 
bly suspect. Sources #3 & 4 were not 
found in the catalog, consistent with the 
lowest significance as reported in Table 2 
and below the significance threshold of the 
ROSAT 1RXH catalog. As an additional 
check, we used the CIAO tool wavdetect 
to search for sources in the ROSAT HRI 
image. Only one source was found, cor- 
responding to source #2, with a signifi- 
cance of 4.1 a 11 consistent with detecting 
just the brightest of the 5 sources. Un- 
fortunately, CIAO does not offer the more 
sensitive PSF-fitting algorithm. 

Of the 5 ROSAT sources in Table 2, only 
CXOU J205847.5+414637 (source #2) was 
detected with Chandra. Source #1 fell 
between the ACIS-I and ACIS-S chips. 
Sources #3-5 were marginal ROSAT de- 
tections and unless they were very soft or 
transient, should have been detectable by 
Chandra. Hence they are likely spurious. 

If the energy spectrum is assumed the 
same for the Chandra and ROSAT obser- 
vations of CXOU J205847.5+414637, the 

The significance was found by dividing the 
net counts by the "Gehrels error" of the back- 
ground counts, cr = l + (BKG.COUNTS + 
0.75) 1/2 . See the Detect Manual on 

http: / / asc.harvard.edu/ ciao/manuals.html for 
details. 
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Table 2 

Sources Detected with ROSAT 



# 


R.A. 


Dccl. 


Total Counts Rate 


Error 


Catalog Sources a 








(cts s" 1 ) 


(cts s" 1 ) 





1 314.7504° 41.8786° 11.1 0.0049 0.0016 1RXH J205900.2+415242 

1RXH J205900.0+415242 

2 314.6972° 41.7771° 13.5 0.0058 0.0016 1RXH J205847.5+414637 

1RXH J205847.3+414638 
1RXH J205847.3+414638 
1RXH J205847.5+414641 

3 314.9085° 41.7064° 6.3 0.0027 0.0012 

4 314.7629° 41.6315° 4.3 0.0019 0.0010 

5 314.8119° 41.6086° 7.3 0.0032 0.0013 1RXH J205915. 0+413631 



a Sources within errors from "The First ROSAT Catalog of Pointed Observations with the High 
Resolution Imager (ROSATHRITOTAL / 1RXH)" available through http://heasarc.gsfc.nasa.gov 
and http: / / www.xray.mpe. mpg.de/cgi-bin/rosat/src- browser. 



ROSAT count rate corresponds to about 
10 times brighter than that observed with 
Chandra, indicating that the flux varied 
significantly between the two observations. 
This is similar to GRO J2058+42 (see Fig- 
ure 1) which was actively outbursting in 
1997, while the outbursts faded below de- 
tectability with RXTE by the time of the 
Chandra observation. 

4. CXOU J205847.5+414637 Opti- 
cal/IR Observations and Results 

4.1. Archival Searches 

We used BROWSE 12 to search for cat- 
aloged objects within the 99% confidence 
the X-ray error circles listed in Table 1. 
Counterparts were found only for CXOU 
J205847.5+414637 in the USNO-B1.0 and 



see http : / /heasarc . gsfc . nasa.gov/db- 

perl/W3Browse/w3browse.pl 



2MASS catalogs. The angular separa- 
tion between these two objects is 0.13". 
The 2MASS object magnitudes were: J = 
11.740 ± 0.022, H = 11.282 ± 0.018, and 
K = 10.930 ± 0.017. The 4' RXTE error 
circle contained 605 USNO-B1.0 sources 
and 636 2MASS sources. The correspond- 
ing densities of sources were used to calcu- 
late the expected average number of acci- 
dental coincidences, A^ 99 , listed in column 
5 of Table 3. The probability of getting 
one or more matches by chance is given by 
the Poisson probability, 1 — e~ Nr ", which 
for small values of the exponent is approx- 
imately N r99 . 

4.2. Optical Photometric Observa- 
tions 

The field around the best-fit RXTE po- 
sition for GRO J2058+42 was observed 
through the B, V, R, and I filters and a 

narrow filter centered at 6563 A (Ha filter) 
using the 1.3 m telescope of the Skinakas 
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Table 3 

Candidate Counterparts to the Chandra X-ray Sources 



Chandra Name 


Catalog 


R.A. 
(J2000) 


Dec. 
(J2000) 


N r99 * 


101 


USNO 






0.19 




2MASS 






0.20 


CXOU J205847.5+414637 


USNO 


314.698078° 


41.777028° 


0.021 




2MASS 


314.698057° 


41.777000° 


0.022 


b 




314.6979° 


41.7769° 


0.001 



a The average number of accidental coincidences expected in the 99% confi- 
dence Chandra error circle. 

b From our optical observations. 



observatory on several occasions through- 
out summer 2003 and 2004 (see Table 4). 
The telescope was equipped with a 1024 x 
1024 SITe CCD chip with a 24 ^m pixel 
size (corresponding to 0.5" on the sky). 
The telescope was pointed at the RXTE 
position R.A = 314.75°, Decl. = +41.72° 
and the field of view was ~ 8.4' x ~ 8.4'. 
Standard stars from the Landolt (1992) 

and Oja (1996) lists were used for the 
transformation equations. Reduction of 
the data was carried out in the standard 
way using the IRAF tools for aperture pho- 
tometry. The results of the photometric 
analysis are given in Table 4. 

Figure 3 shows a Johnson i?-band im- 
age of the field around GRO J2058+42 ob- 
tained from the Skinakas 1.3-m telescope 
on 2003 June 8. The RXTE error cir- 
cle and the optical counterpart to CXOU 
J205847.5+414637 have been marked. The 
second Chandra source, 101 in Table 1, and 
the lower portion of the RXTE error circle 
were outside the field of view. 

An optical color-color diagram was gen- 
erated with the instrumental magnitudes 
of the B, V, R, and Ha filters obtained 
during the 2003 June 8 and 2004 July 5 
observations by plotting the "red color" 
R — Ha as a function of the "blue color" 




Fig. 3. — 8.4' x 8.4' R-band image taken 
from the Skinakas Observatory on 2003 
June 8. The RXTE error circle, centered 
on a = 314.75°, 5 = +41.72° (J2000), 
and the proposed optical counterpart, a = 
314.6979°, 5 = +41.7769° (J2000), are 
shown. Our 99% confidence 1.4" radius 
Chandra error circle, centered at a = 
314.69794°, 5 = +41.77704° (J2000) falls 
entirely on the marked optical star in this 
image. Source 101 and the bottom edge 
of the RXTE error circle were outside the 
field of view. 

B — V (See Figure 4). Be stars showing 
emission in Ha are expected to occupy the 
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Table 4 

Photometric Measurements of CXOU J205847.5+414637's Counterpart 



Date 


MJD 


B 


V 


R 


I 


2003 June 7 


52798.50 


16.01+0.03 


14.91±0.03 


14.22+0.03 




2003 June 8 


52799.44 


16.07+0.03 


14.94±0.02 


14.24+0.02 




2003 August 24 


52876.37 


16.03+0.02 


14.90±0.03 


14.25+0.03 


13.49±0.03 


2004 July 5 


53192.46 


16.03+0.02 


14.88±0.02 


14.16+0.02 


13.35±0.02 


2004 July 27 


53214.51 


16.07+0.02 


14.92±0.03 


14.19+0.03 


13.41±0.04 


2004 August 24 


53242.43 


16.05+0.02 


14.95±0.02 


14.24+0.02 


13.46±0.02 



upper left region of the diagram. Be stars 
are expected to show low (blue-dominated) 
B — V colors because they are early-type 
stars (although they normally appear red- 
der than non-emitting B stars due to the 
circumstellar disk) and also bright (less- 
negative) R — Ha colors because they are 
relatively strong Ha emitters. 

The color-color diagram taken alone 
provides two potential candidates to be 
the optical counterpart associated with 
GRO J2058+52. The two stars lie out- 
side the 90% RXTE error circle, but while 
the data point marked with a star-like 
symbol (a = 314.6979°, 5 = +41.7769°, 
J2000) is 18" outside the error circle and 
corresponds to the X-ray source CXOU 
J205847.5+414637 detected with Chandra, 
the other candidate (the one showing the 
largest R — Ha color values and located at 
a = 314.6833°, 5 = +41.7878°, J2000) is 
about 1.2' away from the error circle and 
was not detected with Chandra or ROSAT. 
Also, the optical spectrum of the other 
candidate does not resemble that of an 
early-type star. 



in 




1 .5 2 

B-V 



Fig. 4. — Optical color-color diagram of 
the field around GRO J2058+42. The data 
point marked with a star symbol is the 
CXOU J205847.5+414637 optical counter- 
part. 

4.3. Optical Spectroscopic Observa- 
tions 

Optical spectroscopic observations of 
the CXOU J205847.5+414637 counterpart 
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were obtained from the Skinakas obser- 
vatory in Crete (Greece) on 2004 June 
25, July 6, and August 26, and from the 
William Herschel Telescope (WHT) lo- 
cated at the Roque de Los Muchachos ob- 
servatory in La Palma (Spain) on 2004 July 
4. The 1.3 m telescope of the Skinakas Ob- 
servatory was equipped with a 1024x1024 
Thomson CCD and a 1302 1 mm -1 grat- 
ing, giving a nominal dispersion of ~1.3 
A/pixel. The WHT instrumental set-up 
during the service run utilized the ISIS 
double-arm spectrograph. The blue arm 
used the R1200B grating with wavelength 
coverage 3890 - 4440 A, while the red arm 
used the R1200R grating with a wave- 
length coverage of 6180 - 6920 A. The spec- 
tra were reduced using the STARLINK 
Figaro package (Shortridge et al. 2001) 
and analyzed using the STARLINK Dipso 
package (Howarth et al. 1998). 

These spectroscopic observations re- 
vealed Ha emission with split profiles. 
Equivalent widths (EW) for the Ha pro- 
files are listed in Table 5. In the red band 
WHT spectrum, top panel in Figure 5, 
the He 6678 A line also shows a double- 
peaked profile. The split profiles provide 
further evidence for a Be star, since the 
double-peaked shape can be interpreted as 
coming from the circumstellar envelope. 
The blue band spectrum, lower panel in 
Figure 5 showed various other members of 
the Balmer series in absorption. 

The blue band spectrum lacks, however, 
strong He I lines, such as 4009, 4144, and 
4387 A. Only He I 4026 and 4121 A (this 
one was rather weak) are present. The 
strength of He I lines reaches a maximum 
for spectral type B2 and falls off on either 
side (i.e., Bl, B3, etc). Thus, we are likely 



- DIB DIB DIE 






^800 4000 4200 4400 

Wavelength (A) 

Fig. 5. — The red (top) and blue (lower) 
spectrum of the counterpart to CXOU 
J205847.5+414637 taken on 4 July 2004. 
A strong, split Ha profile is visible with 
an equivalent width of —4.5 ± l.oA. The 
He 6678A line is also clearly visible, 
various Diffuse Interstellar Bands. In the 
blue, several lines from the Balmer series, 
the He I lines at 4026A and 412lA and Si 
IV 4089A can be seen in absorption. 

dealing with a very early B or late O star. 
The Si IV 4089 A line is quite prominent. 
In fact, the ratio Si IV 4089 A/He I 4121 A 
> 1. This could be explained assuming an 
earlier spectral type (earlier than 09) or 
else with a more evolved companion. The 
first possibility is ruled out by the lack of 
He II (for example, He II 4200 A). The 
lack of information above 4400 A makes it 
difficult to perform a good spectral classifi- 
cation, but the data suggest a classification 
of O9.5-B0IV-V 
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Table 5 

Equivalent Widths of the Ha Line 



Date 


MJD 


Equivalent Width (A) 


2004 June 25 


53181 


-4.2 ±0.6 


2004 July 04 


53190 


-4.5 ± 1.0 a 


2004 July 06 


53192 


-5.0 ±0.9 


2004 August 26 


53244 


-4.7 ±0.5 



Affected by a cosmic ray hit. 



5. Discussion 

5.1. Distance Estimate from Opti- 
cal/IR data 

In order to estimate the distance it is 
important to make a good estimate of the 
amount of absorption to the source. 

The interstellar absorption to the source 
can be estimated from the strength of 
the diffuse interstellar bands (DIB, Herbig 
1975). The measurement of the equivalent 
widths of the DIB is, however, hampered 
by the low signal-to-noise ratio of the op- 
tical continuum. Using the strongest lines 
(those at 6202 A and 6613 A), the esti- 
mated color excess is E(B — V) = 1.3±0.1, 
but the number of measurements is only 3. 
If we add the 5778/80, 6195 and 6269 A 
lines then the average value of the redden- 
ing is E(B - V) = 1.2 ± 0.2, where the 
error is the standard deviation of 7 mea- 
surements. The main uncertainty in the 
derivation of the equivalent widths comes 
from the difficulty in defining the contin- 
uum. 

This color excess is consistent with that 
estimated from the photometric data. A 



B0V star has an intrinsic color (B — V) = 
—0.26 (Wegner 1994) and taking the mea- 
sured photometric color (B — V) = 1.12 ± 
0.04 we derive an excess E(B — V) — 
1.38 ± 0.04. The slightly higher reddening 
obtained from the photometric magnitudes 
may be explained by the contribution of 
the circumstellar disk around the Be star. 
In contrast, the interstellar lines should be 
free of such effects. 

Using the Ha EW ~ 4.5 A the expected 
IR colors may be estimated for a Be star 
plus circumstellar disk. Using Figure 8 
from Coe et al. (1994) the intrinsic (J — K) 
is shown to be in the range —0.03 ± 0.05. 
This gives E(J - K) ~ 0.84 ± 0.05 and 
a corresponding E(B — V) ~ 1.6 ± 0.1. 
Though slightly higher than the value of 
1.38 ± 0.04 determined from the more pre- 
cise photometry, the agreement is suffi- 
ciently good to confirm that both the op- 
tical and IR data are likely to be from the 
same object. 

Taking into account the above three es- 
timates of the reddening we determine the 
weighted mean value to be E(B — V) — 
1.4 ±0.1. 
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Figure 6 shows the photometric values 
(optical and infrared) de-reddened by E(B- 
V)=1.4. Superimposed on the these val- 
ues is a stellar atmosphere for a BO star 
(T eff = 28,000K and hgg = 4.0) (Kurucz 
1979). The model has been normalized to 
the B band data since it is assumed that 
this band will not be significantly affected 
by any contribution from the circumstellar 
disk. It is clear from this figure that the 
IR fluxes lie significantly above the model 
fit, indicating a strong IR excess. This is 
to be expected from a Be star, though the 
strength of the excess in this case is large 
- perhaps partially because the data were 
not taken contemporaneously. Nonethe- 
less, the general quality of fit to the optical 
fluxes is further confirmation that our esti- 
mate for the reddening is reasonably good. 




1000 1300 1400 1600 1600 

Wavelength (nm) 



Fig. 6. — Fit of Kurucz atmospheric model 
for a B0/B1 star to the optical and infrared 
fluxes de. Squares denote the photometric 
values (optical and infrared) de-reddened 
by E(B — V) = 1.4. The solid line denotes 
a stellar atmosphere model for a BO star. 



Taking the standard law Ay = 3AE(B— 
V) Rieke & Lebofsky (1985) and as- 
suming an average absolute magnitude 
of My = -4.2 ±0.1 (Vacca, Garmany, 
& Shull 1996) the distance to GRO 
J2058+42 is estimated to be ~ 9.0 ± 1.3 
kpc. The error was obtained by propa- 
gating the errors in my, E(B — V) and 
My. Note that the error in the absolute 
magnitude accounts for the uncertainty in 
the spectral type of the optical companion 
(O9.5-B0.5). 

The distance estimate from the opti- 
cal/IR counterpart to CXOU J205847.5+414637 
is consistent with X-ray distance estimate 
for GRO J2058+42 of 7-16 kpc (Wilson et 
al. 1998), based on measurements of spin- 
up and flux in the giant outburst in 1995. 

5.2. Comparison with Other Quies- 
cent Be/X-ray Binaries 

From our Chandra observations, CXOU 
J205847.5+414637's energy spectrum was 
approximately fitted with a power-law 
with a photon index ~ 1.8 and an un- 
absorbed flux -Fi_ 10 kev = (3 — 9) x 10 
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33 



- erg cm 2 s 1 , or Li_io kev = (3 — 9) x 10' 
erg s -1 for a distance of 9 kpc. Campana 



et al. (2002) measured power-law indices of 
2.1 and 2.6 and unabsorbed L .5_io kev — 
(l-3)xl0 35 and (0.8-2) x 10 33 erg s' 1 for 
A0538-66 and 4U 0115+63, respectively, 
with BeppoSAX. No pulsations were de- 
tected for either source, with a 3a upper 
limit of 30% for 4U 0115+63. Chandra de- 
tected 22 photons from V0332+53 in quies- 
cence corresponding to Lo.5-10 kev ~ 10 33 
erg cm -2 s _1 (Campana et al. 2002). No 
optical observations were reported for 
these objects, so the state of the Be disk 
was unknown. Quiescent observations 
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of A0535+26 with BeppoSAX (Orlandini 
et al. 2004) and RXTE (Negueruela et 
al. 2000) are most intriguing. In 1998, 
A0535+26 was detected at L 3 _ 2 o kcv =~ 
4 x 10 33 erg s _1 with RXTE, with a power- 
law index of ~ 3. Pulsations were detected 
during this observation with a pulse frac- 
tion of > 53%. Optical observations of 
V725 Tau, the counterpart to A0535+26, 
showed Ha in absorption from the underly- 
ing star, indicating that the circumstellar 
disk was gone (Negueruela et al. 2000). 
BeppoSAX observations in 2000-2001 also 
revealed pulsations with a ~ 50% frac- 
tion, a power-law index of ~ 2 and an 
unabsorbed L 2 -io kcv = (1-5 - 4.4) x 10 33 
erg s -1 . However optical observations in- 
dicated that Ha emission, i.e., the cir- 
cumstellar disk, had returned. CXOU 
J205847.5+414637's luminosity and spec- 
tral shape are similar to that observed from 
A0535+26, especially in the BeppoSAX ob- 
servations where the Be disk had returned. 
If CXOU J205847.5+414637 and GRO 
J2058+42 are the same object, it is also 
a long-period pulsar similar to A0535+26. 

The RXTE ASM light curve showed 
that regular outbursts of GRO J2058+42 
continued from its discovery in 1995 until 
at least mid-2002. This long period of ac- 
tivity from 1995-2002 was similar to EXO 
2030+375 (Wilson et al. 2002). For EXO 
2030+375, the extended activity was be- 
lieved to be the result of the Be disk being 
truncated at a particular resonance radius 
(4:1), so that an outburst was expected 
at every periastron passage unless the Be 
disk disappeared (Okazaki & Negueruela 
2001). After mid-2002, GRO J2058+42's 
outbursts faded below detectability in the 
RXTE ASM. If CXOU J205847.5+414637 



and GRO J2058+42 are the same ob- 
ject, our optical observations show that 
although GRO J2058+42's outbursts ap- 
peared to have ceased by 2004, detec- 
tion of Ha emission indicated that the 
Be disk was still present. Hence GRO 
J2058+42's behavior, while similar to EXO 
2030+375's during outbursts, does not fit 
the EXO 2030+375 model predictions for 
outburst cessation. EXO 2030+375's out- 
bursts have not ceased to date, so it is not 
clear whether or not EXO 2030+375 will 
fit the model either. 

6. Conclusions 

Clearly we have discovered a Be/X-ray 
binary CXOU J205847.5+414637. Differ- 
ent estimates of the number of Be stars 
in the Galaxy give radically different num- 
bers, but all give a low chance probability 
of finding a serendipitous Be star in our 
Chandra error circle. Based on evolution- 
ary arguments, Meurs & van den Heuvel 
(1989) predict 2000-20,000 Be/X-ray bina- 
ries in the galaxy. Porter & Rivinius (2003) 
state that about 1/3 of known Be stars are 
in binaries. About 2/3 of these binaries 
are Be/X-ray binaries, i.e., they contain a 
neutron star. This leads to roughly 10,000- 
100,000 Be stars in the Galaxy. If we as- 
sume an angular size of 360° x 20° for the 
Galactic plane, that implies 0.004-0.004 Be 
stars per arcmin 2 . Hence, the chance prob- 
ability of finding an unrelated Be star in 
our Chandra 99% confidence error circle is 
0.00007-0.0007%. Alternatively, the mass 
of the Galaxy is 1.8 x 10 n M o (Zombeck 
1990), meaning about 10 11 stars in the 
Galaxy. About 1 star in 800 is a B star 
(Binney & Merrifield 1998) and about 
17% of B stars are Be stars (Zorec & Briot 



15 



1997). This gives us a much larger number 
of Be stars in the Galaxy, about 2 x 10 7 , 
leading to about 0.8 Be stars per arcmin 2 
and a chance probability of 0.1% of find- 
ing an unrelated Be star in our Chandra 
error circle. Neither argument takes into 
account the fact that not all Be stars show 
Ha emission at a given time and some per- 
centage of them will not be observable due 
to absorption effects, so the actual proba- 
bility is even lower. In addition, Be stars 
in Be/X-ray binaries cover a very narrow 
range in spectral type (09-B2, Negueru- 
ela 1998) and the spectral type of CXOU 
J205847.5+414637's companion is in that 
range. 

Based on evolutionary arguments (Meurs 
& van den Heuvel 1989), we should expect 
to find 0.004-0.04 Be/X-ray binaries in the 
RXTE error circle for GRO J2058+42 and 
0.02-0.2 in the ACIS-I field. However, 
based on our alternative argument, we 
would expect to find about 8 Be/X-ray bi- 
naries in the RXTE error circle and about 
44 in the ACIS-I field. Both estimates 
assume that all Be/X-ray binaries in the 
galaxy are detectable. To our knowledge 
to date no previously unknown quiescent 
Be/X-ray binaries have been discovered in 
X-rays, suggesting that either they are less 
common than predicted or are too faint 
to observe. Our observations and observa- 
tions of other Chandra fields (e.g., Rogel 
et al. 2004), suggest that the alternative 
argument overestimates the number of Be 
stars. Unfortunately little is known about 
quiescent Be/X-ray binaries. 

Because pulsations were not detected 
with our Chandra observation, we cannot 
definitively say that CXOU J205847.5+414637 
is GRO J2058+42. However, beyond the 



positional coincidence, other evidence is 
suggestive that they are the same source. 
CXOU J205847.5+414637 was the bright- 
est object observed with ROSAT in 1997, 
when GRO J2058+42 was active, and 
it was roughly 10 times brighter than 
in the Chandra observations. Similarly, 
RXTE PCA observations showed that 
GRO J2058+42 was a factor of about 
10-100 brighter in outburst in 1998 than 
upper limits in 2003 December. The 
ROSAT observations approximately cor- 
responded in orbital phase to the faintest 
RXTE detections. Longer-term obser- 
vations with the ASM also showed that 
GRO J2058+42 had faded. Further, GRO 
J2058+42 shows classic Be/X-ray binary 
behavior and CXOU J205847.5+414637 is 
associated with a Be star. Lastly, the 7- 
16 kpc distance to GRO J2058+42 is in 
agreement with the distance of 7.7-10.3 
kpc estimated for the optical counterpart 
to CXOU J205847.5+414637. Additional 
observations of CXOU J205847.5+414637 
are needed to be certain if it is a new qui- 
escent Be/X-ray binary or if it is GRO 
J2058+42. 

This publication makes use of data 
products from the Two Micron All Sky 
Survey, which is a joint project of the Uni- 
versity of Massachusetts and the Infrared 
Processing and Analysis Center /California 
Institute of Technology, funded by the 
National Aeronautics and Space Admin- 
istration and the National Science Foun- 
dation. In addition, this research has 
made use of data obtained from the High 
Energy Astrophysics Science Archive Re- 
search Center (HEASARC), provided by 
NASA's Goddard Space Flight Center 
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(GSFC) and quick-look results provided 
by the RXTE/ ASM teams at MIT and at 
the SOF and GOF at GSFC. 
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